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ABSTRACT 

The discovery that a substantial fraction of Kuiper Belt objects (KBOs) ex- 
ists in binaries with wide separations and roughly e qual masses, ha s moti vated 



a variety of new theories explaining their formation. iGoldreich et al.l (120021 ) pro- 
posed two formation scenarios: In the first, a transient binary is formed, which 
becomes bound with the aid of dynamical friction from the sea of small bodies 
{LP'S mechanism); in the second, a binary is formed by three body gravitational 
deflection (L^ mechanism). Here, we accurately calculate the L^s and Lr" forma- 
tion rates for sub-Hill velocities. While the Li^s formation rate is close to previous 
order of magnitude estimates, the formation rate is about a factor of 4 smaller. 
For sub-Hill KBO velocities {v <^ vh) the ratio of the to the L^s formation 
rate is 0.05(v/vjy) independent of the small bodies' velocity dispersion, their sur- 
face density or their mutual collisions. For Super-Hill velocities {v ^ vh) the 

mechanism dominates over the LP's mechanism. Binary formation via the LP 
mechanism competes with binary destruction by passing bodies. Given sufficient 
time, a statistical equilibrium abundance of binaries forms. 

We show that the frequency of long-lived transient binaries drops exponen- 
tially with the system's lifetime and that such transient bin aries are n o t imp ortant 



for binary formation via the LP mechanism, contrary to iLee et al.l (120071 ) . For 
the LPs mechanism we find that the typical time, transient binaries must last, to 
form Kuiper Belt binaries (KBBs) for a given strength of dynamical friction, Z), 
increases only logarithmically with D. Longevity of tr ansient binaries (with life- 



times ^ 151] ^ as suggested by lAstakhov et al.l (120051 )) only becomes important 



for very weak dynamical friction (i.e. D < 0.002) and is most likely not crucial 
for KBB formation. 



Subject headings: Kuiper Belt — planets and satellites: formation 
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1. INTRODUCTION 

One of the many intriguing discoveries in the Kuiper Belt is that a substantial fraction of 
its largest members a r e bina ries. 48 such systems are currently known (for a comprehensive 



review see 



Noll et al.l (l2007l )). Broadly speaking, we can identify two classes of Kuiper Belt 
binaries (KBBs). The first class consists of small satellites around the largest Kuiper Belt 
objects (KBOs) and the second of roughly equal-mass binaries with wide separations. The 
existence of the first class of binaries is most likely explained by the standard formation 
scenario involving a collision and tidal ev olution, as has been proposed for the formation of 



the Moon and th e Pluto-Charon system (IHartmann fc David Il975l : ICameron fc Wardlll976 



McKinnonI Il989l ) . This formation scenario fails however for the second class of KBBs, since 



the formation of comparable mass KBBs (e.g. 


Weidenschi 


ling 


2002; 


Gold 


reich et al. 


2002; 


Funato et al. 


2004; 


Astakhov et al. 


2005; 


Lee et al. 


2007 


)• 


Weidenschilline 


( 


2OO2I) proposed 



a new formation mechanism for KBBs consisting of a collision between two bodies inside 
the Hill sphere of a third. However, in the Kuiper Belt, gravitational scattering between 
the two intruders is about 100 time^j more common than a collisio n. Binary forma t ion b y 



three body gravitational deflection (L^ mechanism), as proposed bv lGoldreich et al.l (120021) . 
should therefore dominate over such a collisional formation scenario. iGoldreich et al.l (120021 ) 
proposed a second binary formation scenario: it consists of the formation of a transient 
binary, which becomes bound with the aid of dynamical f rictio n from the sea of small bodies. 



This is called the L'^s mechanism. I Astakhov et al.l (120051 ) and lLee et al.l (120071 ) suggest that 



transient binaries that spend a long time in their mutual Hill sphere, near a periodic orbit, 
form the binaries in the L'^s and mechanism. We address and investigate the relative 
importance of these long-lived transient binaries for the L'^s and Li^ formation mechanism 
and find that they ar e most likely not si gnificant for the overall binary formation in the 



Kuiper Belt. Finally, iFunato et al.l (120041 ) proposed a binary formation mechanism which 
involves a collision between two large KBOs which creates a small moon. An exchange 
reaction replaces the moon with a massive body with high eccentricity and large semi-major 
axis. 

In this paper, we accurately calculate the L'^s and formation rates for sub-Hill KBO 
velocities and discuss how these rates are modified for super-Hill velocities. This allows us 
to determine for which physical parameters and velocity regime each mechanism dominates 
the binary formation. Further, we calculate the frequency of long-lived transient binaries 



^For this estimate we used a 10 ^ and assumed that the velocity dispersion of the KBOs at the time 
of binary formation is less than their Hill velocity, see §2 for details 
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and assess their importance for the overall KBB formation. 

Our paper is structured as follows: In §2 we outline our assumptions, explain our choice 
of parameters and define variables that will be used throughout this paper. We calculate the 
Lr" and L'^s formation rates for sub-Hill KBO velocities in §3 and §4 respectively. We compare 
the L^s and formation rates in the sub-Hill velocity regime in §5. In §6 we discuss how 
these formation rates are modified for super-Hill KBO velocities. The frequency of long-lived 
transient binaries and their significance for the overall KBB formation is calculated in §7. 
Summary and conclusions follow in §8. 



2. DEFINITIONS AND ASSUMPTIONS 



The Hill radius denotes the distance from a KBO at which the tidal forces due to the Sun 
and the gravitational force due to the KBO, both acting on a test particle, are in equilibrium. 
It is given by 

where a is the semi-major axis and M the mass of the KBO . M^, is the mass of the sun. 
We use the 'two-group approximation' (iGoldreich et al.l |2002| . |200J) which consists of the 
identification of two groups of objects, small ones, that contain most of the total mass with 
surface mass density a, and large ones, that contain only a small fraction of the total mass 
with surface mass density S ^ a. We assume a ~ 0.3g cm^^ which is the extrapolation of 
the minimum-mass s olar nebular to a helioce r itric distance of 40A U. Estimates from current 
Kuiper Belt surveys ( iTrujillo fc Brown! l2003l : iTrujillo et al.ll200l[ ) yield E ~ 3 x lO^^g cm~^ 



for KBOs with radii of i? ~ 100 km. We use this value of S, assuming that S during the 
formation of KBBs was the same as it is now. Our choice for E and a is also consistent with 



results from numerical coagulation simulations by lKenyon fc Luul (119991 ). 



Large bodies grow by the accretion of small bodies. Large KBOs viscously stir the small 
bodies, increasing the small bodies' velocity dispersion u. As a result u grows on the same 
timescale as R provided that mutual collisions among the small bodies are not yet important. 
In this case, u is given by 



u 



Vh 



1/2 



(2) 



where a = R/Rh ~ 10""^ at 40AU (IGoldreich et al.ll2002f ). vh is the Hill velocity of the 
large bodies which is given by vh = ^Rh where fl is the orbital frequency around the sun. 
The velocity v of large KBOs increases due to mutual viscous stirring, but is damped by 
dynamical friction from the sea of small bodies such that v < u. Balancing the stirring and 
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damping rates of v and substituting for u from equation [2], we find 



V _2 / S 



— ~ a 



0.1. 



(3) 



For our choice of parameters, we have that v < vh during the epoch of formation of bodies 
with R ~ 100km. Additionally, we argue that v could not have exceeded vh significantly 
during satellite formation in the Kuiper Belt: If v^sc > v > vh, where Wesc is the escape 
velocity from the large bodies, then the timescale for mutual collisions is 



Tcoll ~ 0.13 



3 X 10-4g cm- 



P 



Ig cm 



V 100km J 



a 



1 X 10-4/ yvH, 
1 

Gyr. 



7.9 X 10-iOs-i 



(4) 



Equation H] shows that the collision timescale among the largest KBOs {R > 100km) would 
have been excessively long if v ^ vh during satellite form ation. The ubiquity of small satel- 
lites around KBOs, that have radii as large as ~ 1000km , (IBrown et al.ll2006l : iBrown fc Suer 



20071 ) and the Pluto-Charon system (IWeaver et al.ll2006l ) suggests that v < vh during their 



formation, since t heir origin is best explained by a giant impact (e.g. IStern et al.l 12006 



Brown et al.l 120071 ). This is supported further by the recent discovery of a collisional family 



belonging to ELgi (IBrown et al. 



20071 ). We therefore focus our work on the shear-dominated 



velocity regime {v < vh)- However, we discuss how our results would be modified iiv > vh- 



3. FORMATION RATE 



A transient binary forms when two large KBOs penetrate each other's Hill sphere. This 
transient binary must lose energy in order to become gravitationally bound. In the 
mechanism the excess energy is carried away by an encounter with a third massive body. We 
calculate the binary formation rate via the mechanism in the shear-dominated velocity 
regime. Since the growth of inclinations is suppres sed in the shear-dominated velocity regime 



the disk of KBOs is effectively two-dimensional (IWetherill &: StewartI Il993l : iRafikovl 12003 



Goldreich et al.ll2004l ). We therefore restrict this calculation to two dimensions. As initial 
condition, we assume that all bodies are on circular orbits. We chose to work in the rotating 
frame with the x-axis pointing radially outwards and the y-axis in the prograde direction. 
For a gravitational deflection of three equal-mass bodies, the formation rate per body is 



FRvi 



oo J b2>bi J b\=0 



^pR^ 



-biQFL3{bi,b2,'y)dbidb2d'y. 



(5) 
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S/(47rpi?^/3) is the surface number density of the KBOs. bi and 62 are the relative initial 
separations in the x-direction between bodies 1 and 2 and bodies 1 and 3 respectively. 7 is 
the offset in the y-direction body 3 would have when bodies 1 and 2 would encounter each 
other had their relative velocity been solely due to the Kepler shear of the disk: 3biQ/2. 
Finally, ^^3(61, 62, 7) is a function that takes on the value 1 if the encounter resulted in the 
formation of a binary between any two of the three KBOs involved and otherwise. The 
choice of limits on the integrals in equation [5] ensures no double counting of the binaries. 
Expression can be written as 



FR 



L3 



pR 



where 



27 

32^ 



b2>bi 



00 / ^ 

i^L3(&l,&2,7) 

61=0 



dhi 
Rh 



dh2 
Rh 



d'-f 
Rh 



(6) 



(7) 



Expression ([6]) agrees with the order of magnitude estimate of iGoldreich et al.l (120021 ) if we 
set Ai3 = 1. It is the value of the constant Aj^a we determine here. Since we are interested 
i n close encounters among the KBOs, their interaction is w ell described by Hill's equations 
(jHilllll878l : IGoldreich fc Tremainelll980l : iPetit fc Henonlll986l ) that we modify to include three 
equal mass bodies besides the Sun. The equations of motion, with length scaled by Rh and 
time by for body 1 are given by 



xi — 2yi — 3xi 



jji + 2xi 



3(xi - X2) 



3(xi - X3) 



{{xi - X2)^ + {yi 
3(1/1 - Z/2) 



{{xi - x^y + {yi - y^Yy/'^ 

3(yi - 2/3) 



((xi - X2f + {yi - 2/2)^)^/^ {{xi - X3)2 + (yi - 1/3)2)3/2 ■ 



(9) 



The subscripts 1, 2 and 3 label the x- and y-coordinates of KBO 1, 2, and 3 respectively. 
Similar equations of motion can be obtained for bodies 2 and 3. F/;,3(6i, 62, 7) is calculated 
by numerically integrating the equations of motion. A binary formation event is detected 
in the following way: The equations of motion of the three bodies are integrated until a 
time that corresponds to a separation of at least 30i?// between all three bodies (after their 
conjunction), assuming that their relative velocity is solely due to their Keplerian sheer (i.e. 
ignoring the actual gravitational interaction between the bodies), plus an additional time 
of 120f2~^. If after this time the separation between two bodies is still less than "iRn a 
binary is considered to have formed. We chose a separation of "iRn instead of Rh to allow 
for binary orbits that reach slightly outside Rh- Numerical integrations are terminated 
early if the separation between KBOs becomes less than 10~^i?ii/ and these events are not 
counted towards the binaries formed. This serves two purposes: first of all, IQ'^^Rh roughly 



corresponds to the separation at which physical colhsions occur in the Kuiper Belt. Secondly, 
by introducing a minimum separation, we prevent divergence in the equations of motion. 
This cut-off limits, strictly speaking, the validity of the value of calculated here to 
binary formation at heliocentric distances of ~ 40 AU since the separation in units of Rh, 
corresponding to collisions among the KBOs, is inversely proportional to the heliocentric 
distance. In order to determine A^a we need to cover the three dimensional parameter space 
spanned by bi, 62 and 7. We chose a spacing of O.IRh for all three parameters. 12.5Rh is 
chosen as the upper limit for bi and 62, the upper limit for I7I is 25Rh- The given limits and 
resolution require numerical integrations of ~ 4 x 10^ orbits. We obtain 

A^a = 0.28 ± 0.01 (10) 

where 0.01 is the estimated Poisson error. We repeated the calculation for Aj^a with randomly 
chosen grid points for bi, 62 and 7 and the same number of numerical integrations and 
confirmed that the value of A^a is insensitive to the grid points chosen. A^a tends to 0.35 
in the limit that the bodies are treated as point masses (i.e. the limit that the cut-off tends 
to zero). We will use A^^a = 0.28 since it corresponds to the physically relevant situation in 
the Kuiper Belt. This yields a binary formation rate of 

Fi?i3 = (6.3 ± 0.2) X 10-« f ^ Y( P) Y— ' 

(11) 



1 X 10-4 / V 7.9 X 10-iOs 



which is smaller by l/A^a ~ 4 than the order of magnitude estimate of iGoldreich et al. 



(I2OO2I). 



4. L'^s FORMATION RATE 

So far, we have only considered binary formation due to an encounter with a third 
body that carries away the excess energy. However, binary formation might also occur due 
to dynamical friction generated by the sea of small bodies (L^s mechanism). The random 
velocity of large KBOs is damped due to gravitational interactions with many small bodies. 
Since it is not feasible to examine the interactions with each small body individually, their 
net effect is modeled by an averaged force which acts to damp the large KBOs' non-circular 
velocity. We parameterize the strength of the damping by a dimensionless quantity D defined 
as the fractional decrease in non-circular velocity due to dynamical friction over a time Q^^: 
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The first expression is simply an estimate of dynamical friction by the sea of small bodies 
assuming u > vh- The second expression describes the mutual excitation among the large 
KBOs for V < Vh- These two expressions can be equated since the stirring among the large 
KBOs is balanced by the damping due to dynamical friction. In fact, if v is defined as 
the product of the median eccentricity and the orbital velocity, we can calculate the exact 
relationship between D and (v/vh) since the velocity distribution i n the shear-do r ninat ed 
velocity regime has been fully determined (see ICoUins fc Saril (120061 ) : ICoUins et al.l (120071 )). 
Defining v as the product of the median eccentricity and the orbital velocity, we obtain 



pR \vh 



-1 



(13) 



For p ~ Ig cm~^ and our estimates for {v/vh), S and R from §2 we find D ~ 0.12. We 
calculate the binary formation rate for equal mass bodies via the L^s mechanism in shear- 
dominated velocity regime. As in §3, we restrict this calculation to two dimensions with 
circular motion as initial conditions for the large KBOs and use the same coordinate system 
as in §3. The binary formation rate per body via the L^s mechanism can be written as 



FRl^s 



6=0 



pRy 2 



-bnFL2,{D,b)db 



(14) 



where E/(47rpi?^/3) is the surface number density of the KBOs and b is the relative initial 
separation in the x-direction between the two KBOs. FL2g[D, b) is a function that takes on 
the value 1 if the encounter resulted in the formation of a binary for a given D and b and 
otherwise. Equation [13] can be written as 



FR 



where 



A 



D- 



9 



in 



A 



-2 



a 



FL2s{D,b) 



6=0 



_b_ 

Rh 



dh_ 

Rh 



(15) 



(16) 



Goldreich et al.l (|2002| ) showed, using numerical integrations, that FRi2g is indeed propor- 
tional to D. Here we want to determine the actual value of Ai2g. In Hill coordinates the 
equations of motion of the two KBOs can be decomposed into their center of mass motion 
and their relative motion with respect to one another. The relative motion of two equal mass 
KBOs, including the dynamical friction term, is governed by 

6x 



X — 2y — ?,x 



2x 



6y 



(X2 + ?/2)3/2 



Dx 

x"^ + y2)3/2 

D{y + 1.5x). 



(17) 
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where x and y correspond to the relative separation between the two KBOs in the x and y 
direction respectively. Again, length has been scaled by Rh and time by Vt^^ . Equations [T71 
and [TS] are integrated for different values of D and impact parameters ranging from 2.2Rh 
to 3.2Rh- Impact parameters outside this range result in a distance of closest approach 
between the two KBOs of more than Rh- Figure [T] shows that the rate of binary formation 
is proportional to D. The value of Ai2s, estimated from the line of best fit, is 1.4. This 
yields a binary formation rate of 



FRl2, = 1.3 X 10" 



0.12 



gem \lg cm \100kmy 

ixio-v V7-9 X lo-ios-iy ■ 

i of lOoldreich et all toQ'^ . Although, 



(19) 



Using equation [12] we can retrieve the scalings of iGoldreich et al.l (120021 ). Although, we know 
the exact L'^s formation rate for a given D and have an exact expression for D in terms of 
V (see equation [13]), the relation between v and the actual physical parameters, i.e. the 
numerical coefficient in equation l3l which is needed for a precise value of D is uncertain to 
a factor of order unity. 



Contrary to claims by ]Astakhov et al.] (12005] ) and ]Lee et al.] (12007] ). the L^s mechanism 
does predict a mass-ratio selection. This can be seen from the first part of equation [12] 



For a given u, we have that D oc R^ since vhC(' 



-1/2 



oc R. Large KBOs experience 



stronger dynamical friction (larger D) th an smaller ones. This is not at all surprising an d 
is a general feature of dynamical friction (jChandrasekhar] ]l943] : ]Binney &: Tremaine]]l987] ). 
We can write D = D^M where Dq ~ aC^/u^. Dq is a constant independent of the KBO 
mass for a given a and u. For two KBOs with masses Mi and M2, the position of body 1 
essentially coincides with the center of mass of the two bodies provided that M1/M2 ^ 1. 
In the limit that the KBOs random velocity tends to zero and that DqM2 ^ 1 we can place 
body 1 at the origin of the Hill coordinate system and treat the center of mass as stationary 
throughout the interaction. In this limit we find that the relative motion of the two KBOs 
is governed by 



X 



2y — 3x 



3x 



{X^ + ?/2)3/2 



2D0M2X 



y + 2x 



31/ 



2DoM2{y + 1.5x) 



(20) 



(21) 



(a;2 + 2/2)3/2 

where length is scaled by Rh of KBO 1 and time is scaled by For extreme mass- 

ratio binaries the relevant strength of the dynamical friction that enters equations ]2U] and 
21] is twice that acting on the small body (i.e. 2DqM2), and significantly less than that 
acting on the large body (i.e. DqMi). The L'^s mechanism therfore favors the formation 
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of comparable mass binaries from the largest available bodies over high mass-ratio ones. 
It is an open question, whether this preference for comparable mass binaries remains after 
the Kuiper Belt mass spectrum during their formation and their survival probability are 
accounted for. 



5. COMPARISON OF L'^s AND FORMATION RATES 

We are now able to compare the binary formation rates for the L^s and mechanism 
for sub-Hill velocities. The ratio of the to L^s formation rates is 

£^ = = 0.05— (22) 

where we substituted for D using the exact relationship from equation [131 It is remarkable 
that this expression depends explicity only on v/vh and is independent of what sets D. It 
is therefore independent of the velocity dispersion of the small bodies, their surface density 
or the importance of collisions among the small bodies. We therefore conclude that for 
V <^ Vh, binaries in the Kuiper Belt formed primarily due to dynamical friction rather than 
three body encounters. Figure [U shows the L'^s and formation rates as a function of D. 
For our estimate of {v/vh) ~ 0.1, we have that FRl-a / FRl2s ~ 0.005. 



6. Super-Hill Velocity: v > vr 

Obviously there is some uncertainty in what the actual values of a and S were during 
binary formation. For a few times larger value of S with a unchanged, we enter the regime 
in which v > (this can be seen from equation [2]). Although it is rather unlikely that 
V ^ Vh during binary formation (see §2) we discuss here briefly how this would affect the 
L'^s and formation rates. 

For V > Vh the velocity dispersion of the large bodies is still set by the balance between 
their mutual stirring and the damping due to dynamical friction generated by the sea of 
small bodies. Therefore, dynamical friction shrinks the orbit of a KBB with a mutual orbital 
velocity 

Vb at a rate 

DQ ~ —a-^n ( — ^ (23) 
pR \vhJ 

where we assume that vb < u. For vb ^ v, binaries are broken up by passing KBOs at a 
rate 



Rbreak sa''^ { — ] — ■ (24) 

pR \vh J \vh/ 
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The ratio of these two rates yields 

DQ (vbV 



break 



(25) 



Since the ratio in equation [25] is < 1 for vb < v, we conclude that KBBs with separations 
Rb > Rh{vh/vY (i.e. KBBs with vb < v) tend to be broken up by passing KBOs. Binaries 
with separations of Rcrit = RnivH/vY or less, tend to survive. The cross section for the 
mechanism is therefore reduced with respect to the sub-Hill case. The probability of having 
a KBO within Rent of a given KBO is {EQ)/{pR^v) R^^.^^ where {TjQ)/{pR^v) is the volume 
number density of KBOs. The flux of KBOs into area R^rit is /{pR^v) f -Rcrit- The 
super-Hill formation rate for tight binaries with separations ~ Rcrit, via the mechanism, 
is therefore 



(see also iNoll et al.l (120071 )). In addition to tight binaries with separations of Rcrit and less, 
there exist a second class of binaries with larger separations. Binaries with separations 
Rb > Rcrit are constantly created and destroyed via the mechanism. KBBs can form 
from two KBOs that approach each other with relative velocity vb ^ v while a third KBO 
removes energy, through gravitational scattering, enabling the KBOs to get bound. Since 
we are selecting bodies with relative velocities ~ vb or less, the number of KBOs that can 
form binaries with separation Rb = Ruivu /vbY is reduced by ~ {vb/vY- The formation 
rate for binaries with separation Rb = Rh{vh IvbY is 



These wider binaries {Rb > Rcrit) have a higher formation rate compared to the tight ones 
which have a separation ~ Rcrit- The ratio of the formation rate (equation 1271) to the 
destruction rate (equation yields an equilibrium abundance of binaries per KBO at any 
given time that is given by 

Nkbb ^_^-2(!!ny(v_H\\ ^28) 



Nkbo pR ^ V ^ \vb 

The number of binaries scales as {Rb/ RhY^'^ ■ Binaries with separation Rb are therefore 
{Rb/ RcritY^'^ ~ {v/vbY times more common than those with separation Rcrit provided 
there is sufficient time for the equilibrium to be established. The same statistical equilibrium 
abundance can be derived using phase space arguments. The phase space number density 
of KBOs is {TjVL) / {pR^v^). The phase space volume corresponding to a binary separation 
Rb and velocity vb is R%v% = R%v%{vh/vb)^- Multiplying the KBO phase space number 
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density by the binary phase space volume yields a statistical equilibrium abundance per 
KBO of 

which is in agreement with the binary abundance derived in (128|1 . Whether any of these 
binaries would survive the dynamical excitation of the Kuiper Belt remains an open question. 

The L'^s mechanism fails in creating binaries with separations ~ Rmt since dynamical 
friction is not able to dissipate sufficient energy for tight binaries to form. Dynamical friction 
is only able to assist in the formation of binaries with wide separations (~ Rh) that form 
from KBOs that happen to approach each other with low relative velocities (~ vh)- This 
reduces the number density of KBOs that can participate in binary formation by a factor of 
~ {vh/vY- In this case, the L^s formation rate is given by 

where we have substituted for D from equation [23] in the last step. These wide binaries face 
the same fate as the wide ones formed via the mechanism in that they will be broken up 
quickly due to scattering of other large bodies. However, the L^s mechanism does not even 
contribute significantly to the binary equilibrium abundance calculated in equation [28] since 
FRl2s(Rb ~ Rh) I FRla{Rb ~ Rh) ~ (vh/v)^ ^ 1- Therfore, the L^s mechanism does not 
play an important role in KBB formation if super-Hill velocities prevail. 

In summary, the mechanism forms tight binaries, that tend to be saved from break 
up, at a rate that is reduced by a factor of (vh/v)^^ compared to the sub-Hill case. In 
addition, the mechanism forms wider binaries {Rb > Rcrit), at a higher rate that is 'only' 
reduced by a factor of {vh/v)^{vh/vb)^ relative to the sub-Hill rate. These wide binaries are 
constantly created and destroyed leading to an equilibrium abundance of binaries that scales 
as (Rb / Rh)^^'^ ■ The L^s mechanism is not important if KBOs have super-Hill velocities. 



7. FREQUENCY OF LONG-LIVED TRANSIENT BINARIES AND THEIR 
SIGNIFICANCE FOR BINARY FORMATION 



Astakhov et al.l (l2005[ ) propose that transient binaries, that spent a time of 15Q^^ 
(~ 600yr at 40AU) or longer in their mutual Hill sphere, near a periodic orbit, are re- 
sponsible for binary formation in the L'^s and mechanism. Here, we determine how the 
frequency of long-lived transient binaries depends on the transient binary lifetime. This 
allows us to quantify the importance of long-lived transient binaries for the overall binary 
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formation. Finally, we address the significance of long-lived transient binaries for the Li^s 
and formation mechanism. 



7.1. Frequency of long-lived Transient Binaries 

First, we assess how common long lived transient binaries are. We integrate equations 
[T7] and [18] without the dynamical friction term and determine the time t^R^ over which the 
separation between the two KBOs is less than "iRn for all KBOs that approach one another 
to Rh and less. We chose to calculate the time the two KBOs spent with a separation of less 
than ?)Rh to allow for orbits that reach slightly outside of Rh but return back to within Rh 
during the encounter. We integrate ~ 10^ orbits in total with impact parameters ranging 
from 2.2Rh to 3.2Rh- Impact parameters outside this range result in a distance of closest 
approach between the two KBOs of more than Rh- As initial conditions, we assume that 
the orbits of the bodies are circular. Figure [2] shows that the frequency of transient binaries 
decreases exponentially with the transient binary lifetime, t^^^ . The line of best fit yields a 
differential transient binary frequency, valid for t^ji^ > IQ'^, of 

3^-^ = 1.0 X 10-(o-25*3«H^)^a-2(]. (31) 

The frequency of transient binaries that spend a time of > 15f2~^ with a separation of less 
than ?)Rh is 3 orders of magnitude smaller than that of short-lived ones with ts/?^ > 3fi~^. 
The analysis discussed here was carried out assuming that the KBOs approach each other 
with relative velocities Vrei < vh- Long-lived transient binaries do not exist for bodies that 
encounter each other at Vrei 3> vh- This can be understood by looking at the Jacobi constant. 
The Jacobi constant in Hill coordinates with length scaled by Rh and time by fl~^ is given 
by 

Cj = 3x' ~z'+ f _ ^2 _ . 2 _ . 2^ ^32) 

KBOs that approach each other with Vrei S> vh at Rh will experience at most one encounter 
before they separate. Evaluation of their Jacobi constant at Rh yields that it is large and 
negative. In order to experience multiple encounters, KBOs must approach each other with 
Vrei ~ Vh at Rh which corresponds to Cj of order unity. Since the Jacobi constant is a 
conserved quantity, we can be sure that no long-lived transient binaries exist for KBOs that 
encounter each other at Vrei ^ vh- Long-lived tr ansient binaries t herefore offer no solution 



to the fine tuning problem, contrary to claims by iLee et al.l (120071 ). For KBOs with a given 
velocity distribution there always exist a few bodies that have Vrei < vh even if v ^ vh- 
Such bodies can give rise to long-lived transient binaries in the same way that they can form 
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wide binaries (see §6 for details), but the frequency of transient binaries due to such bodies 
is reduced by a factor of {vfj/v)^. 



7.2. Importance of long-lived Transient Binaries in the Formation 

Mechanism 



Lee et al.l (120071 ) claim that the probability of binary formation from transient binaries 
with ta/f^ < 2.5f2~^ is extremely small and they therefore include only transient-binaries 
with t^ji^ > 5Q^^ in the main set of their integrations. However, their conclusion, that 
the probability of binary formation from transient binaries with t^^^ < 2.5Q~^ is extremely 
small, is due to a bias in their initial conditions that discriminates against binary formation 
from transient binaries with t^ji^ < 5Q~^. The shortcoming of their analysis is due to the fact 
that they launch the third body from an initia l separatio n > 38 Rh when the first two KBOs 



come 



within a few Rh of each othej^. Since Lee et al. ( 2007 ) select the initial conditions 



for the third body such that it penetrates within 2.5Rh, the largest impact parameter is 
~ 4.5Rh- The minimum time it takes for the third body to come within a few Rh of the 
transient binary is therefore ~ 38i?///(1.5 x 4.5i?//f2) ~ 5.6Q~^. The third body therefore 
only reaches the vicinity of the transient binary for t^ji^ > 5.6Q~^, but it is exactly this 
proximity of the third body th at is requi r ed for binary formation by strong gravitational 



scattering. This explains why iLee et al.l (120071 ) find such a small probability for binary 



formation by transient binaries with tsR^^ < 5.QQ~^. The range of impact parameters that 
lead to binary formation is comparable for short- and long-lived transient binaries. This 
means that the transient binary lifetime is the only advantage long-lived transient binaries 
have compared to short-lived ones, in terms of binary formation likelihood. However, as 
we show in §7.1, the frequency of transient binaries drops exponentially as a function of 
their lifetime. The ratio of the binary formation rate due to short-lived transient binaries 
i^sRH ~ 3i7~^) compared to that due to long-lived ones (t^R^ > 15^2"^), is therefore 

FR(t3R„ > 3fi-i) 3^]-l 10-(0-25x3) 

^ " ~ ~ 200 (33) 

FR{t3R„ > 15^]-l) 15^]-l 10-(0-25xi5) ^'"'"^ 

Although the binary formation rate scales linearly with transient-binary lifetime, the fre- 
quency of transient binaries drops exponentially as a function of its lifetime. Therefore, 
long-lived transient binaries are not important for binary formation via the mechanism. 



^The numerical values stated by 
different definitions of Rh 



Lee et al. 



( 2007 ) are multiplied by a factor of 2^/'^ to compensate for the 
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7.3. Importance of long-lived Transient Binaries in the L'^s Formation 

Mechanism 

In general, KBOs that spend a longer time in the Hill sphere, lose more energy due to 
dynamical friction, and are therefore more likely to be captured. However, they might not be 
responsible for the majority of the binary formation, if the frequency for long-lived transient 
binaries is sufficiently small. To address this question we determine the typical time tTyp-, 
required for a transient binary to become bound with the aid of dynamical friction. We define 
txyp as the time it takes for 50% of all the KBOs, that form a binary, to become bound for a 
given strength of dynamical friction D. txyp is measured from the point at which the relative 
separation between the two KBOs is less than 3Rh- We determine tj^yp in the following 
way. First, we integrate the same equations as in §4 (i.e. equations [17] and [18]). We switch 
off the dynamical friction at different times and continue the evolution of the KBOs until 
t = lOOOfi"^. This process is repeated until we find the time for which dynamical friction 
has to act for 50% of all KBOs, that form a binary, to become bound. A transient binary is 
considered to have become bound when it remains a binary (i.e. relative separation < 3Rh) 
until t = 1000f2~^. We repeat this for different D in order to reveal the relationship between 
txyp and D. Again, impact parameters are chosen to range from 2.2Rh to 3.2Rh. Figure [3] 
shows that, for D > 0.002, the typical time for permanent capture does not depend linearly 
on the strength of the dynamical friction D, but shows a weaker logarithmic dependence. tTyp 
only ranges from ~ 2Q~^ for D ~ 0.2 to ~ 10fi~^ for D ~ 0.002. Furthermore, figure [3] shows 
a noticeable break around D ~ 0.001; for D < 0.001, txyp increase significantly to 20Q~^ 
and more. From this , we conclude that longevity of the transient binary (as discussed by 



Astakhov et al.l (|2005l ) with t^ji^ ^ 15Q^^) becomes only important for very weak dynamical 
friction (i.e. D < 0.002) and is most likely not crucial for KBB formation. In §4 we estimate 
D ~ 0.12, in which case longevity of transient binaries (tsi?^ ^ 15fi~^) is unlikely to be a 
major requirement for binary formation. 



8. SUMMARY AND CONCLUSIONS 

We accurately determine the L^s and formation rates for v < vh- We find that while 
the L^s formation rate is close to previous order of magnitude estimates, the formation 
rate is about a factor of 4 smaller. For v <^ vh-, the ratio of the to the L'^s formation 
rates is ~ 0.05(f /f//) and is independent of what sets D. It is therefore independent of the 
velocity dispersion of the small bodies, their surface density or the importance of collisions 
among the small bodies. For sub-Hill KBO velocities, binaries in the Kuiper Belt formed 
primarily due to dynamical friction rather than three body encounters. For super-Hill KBO 



- 15 - 



velocities {v ^ vh) the LP's mechanism becomes unimportant. The Li^ mechanism forms 
tight binaries that tend to be saved from break up at a rate that is reduced by a factor of 
{vh/vY^ compared to the sub- Hill case. In addition, the mechanism forms wider binaries 
{Rb > Rcrit), at a higher rate that is 'only' reduced by a factor of {vh /v)^{vh /vbY relative 
to the sub-Hill rate. These wide binaries are constantly created and destroyed leading to 
an equilibrium abundance of binaries that scales as {Rb/ RhY^"^ ■ Whether and how any of 
these wide binaries would survive the dynamical excitation of the Kuiper Belt remains an 
open question. 

In addition, we determine the frequency of long-lived transient binaries. We show that 
the frequency of long-lived transient binaries drops exponentially with the system's lifetime 
for Vrei < Vh- About 1000 transient binaries occur with ts/j^ > 30^^ for each transient 
binary with tsB^^ > 150"^. T he long-lived transient binaries investigated by lAstakhov et al. 



(120051 ) and iLee et al.l (120071 ) are therefore very rare. Long-lived transient binaries are not 
important for binary formation via the Li^ mechanism, since the binary formation rate scales 
only linearly with transient-binary lifetime, but the frequency of transient binaries drops 
exponentially as a function of its lifetime. Long-lived transient binaries do not exist for 
Vrei 3> Vh- We show that the apparent shortage of binar ies forming from short-lived transient 
binaries (i.e. t^RH ^ 2.5fi~^) found by iLee et al.l (120071 ) can be explained by a bias in their 
initial conditions that discriminates against binary formation from transient binaries with 
t^Ru ^ 5fi~^. Finally, to assess the importance of long-lived transient binaries in the L^s 
mechanism, we determine the typical time tTyp required for a transient binary to become 
bound with t he aid of dynani i cal fri ction. We show that longevity of the transient binary (as 
discussed by lAstakhov et al.l (120051 ) with t^R^^ ^ 15fi"^) only becomes important for very 
weak dynamical friction (i.e. D < 0.002). We estimate D ~ 0.12, in which case longevity of 
transient binaries {^Rjj ^ 15i7~^) is unlikely to be a major requirement for binary formation. 
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Fig. 1. — Binary formation rate as a function of dynamical friction strength D. The crosses 
correspond to the formation rate via the Li^s mechanism and the dashed horizontal line 
corresponds to the Li^ formation rate for {T,/pR)a~'^ = 3 x 10~^. The L^s formation rate 
is proportional to D. In §4 we estimate that D ~ 0.12 as a result of which FRLs/FRijig ~ 
0.005. 
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Fig. 2. — Differential transient binary frequency A{FRti,) as a function of the transient 
binary lifetime t^R^^ in the shear dominated velocity regime. t^R^ is the time the transient 
binary separation is less than 3Rh- The frequency decreases exponentially with t^Rj^fl. 
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Fig. 3. — The variation of txyp, the time taken for 50% of binaries to get bound, plotted 
against strength of dynamical friction D. For about two orders of magnitude change in D 
{D ~ 0.2 to D ~ 0.002) txyp only changes from ~ 2^"-*^ to ~ 10Q~^. A rapid rise in txyp 
occurs for D < 0.001. 



